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Abstract 

A new approach to summation of divergent field-theoretical series is sug- 
gested. It is based on the Borel transformation combined with a conformal 
mapping and does not imply the knowledge of the exact asymptotic parame- 
ters. The method is tested on functions expanded in their asymptotic power 
series and applied to estimating the ground state energy of simple quantum 
mechanical problems including anisotropic oscillators and caclulating the crit- 
ical exponents for certain comformal field models. It can be expected that the 
new approach to summation may be used to obtaining numerical estimates 
for important physical quantities represented by divergent series in two- and 
three-dimensional field models. 
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In this brief report we discuss the problem of summation of "bad" series arising 
in various fields of physics |], p|. By the strong coupling we mean the situation 
when the parameter of expansion does not belong to the range of convergence of the 
perturbative series. The typical situation is when the radius of convergence is just 
zero and the resulting series are asymptotic. On the other hand, these series can 
represent important physical quantities, and to extract reliable information from 
them, they should be processed by a proper resummation technique. 

The main goal of the present work is to suggest a new approach to treating di- 
vergent series and to apply it to finding numerical estimates of the critical exponents 
for two- and three dimensional anisotropic field models. The resummation technique 
proposed is tested on simple model functions expanded in their asymptotic power 
series, estimating ground state energy of isotropic and cubic anharmonic oscillators, 
Yukawa potential, and critical exponents for conformal field theories. 

At present, there are several resummation procedures such as simple Pade, Pade- 
Borel, and Pade-Borel-Leroy techniques, whose application, however, is limited to 
series with coefficients alternating in signs. For the more sophisticated method based 
on the Borel transformation combined with a conformal mapping, first proposed in 
Refs. P, m, this limitation is not crucial. This technique was then elaborated and 
systematically used for various phase transition problems || §, |7|, §, |9], [II], [11], [L2|, 



13| , |14| being now regarded to as a most universal procedure. But it requires the 
knowledge the exact asymptotic high-order behavior of the series. As a rule the 
coefficients of the series J2kfk9 k behave at large k as k\k b °( — a ) k . The numbers a 
and b Q characterize the main divergent part. Nowadays these parameters are found 
only for the simplest case of the 0(iV)-symmetric models [15], |I6, [17|], and calculating 



them for anisotropic models is a most difficult problem as yet unsolved. As an 
exception we mention the anisotropic quartic quantum oscillator which represents 
a one-dimensional (p 4 field theory with the cubic anisotropy. For the perturbation 
expansion of the ground state energy of this system the asymptotic parameters were 



found in Ref. | lq| . Within the assumption of the weak anisotropy the large-order 



asymptotic behavior of the /^-functions for the cubic model was deduced in Ref. [19 



and then used for determination of the stability of the cubic fixed point in three 



dimensions pC ] 



Below, we make an attempt to overcome the outlined difficulties and suggest a 
new approach to summation of the divergent field-theoretical series, which is based 
on the standard technique Borel transformation combined with a conformal mapping 
Jl], HJ, but which does not involve the exact values of the asymptotic parameters. 
We start from the Borel-Leroy transformation modified with a conformal mapping 
in the form If7fl 



LO k (x) 



The coefficients A^(A) are determined from the equality B(x{u)) = tt^ttsx where 



ag' ^ag' (1 — uj k (x)) 2X 
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fx+X- 



/x+1- 
fk 



j and the Borel transform B(x) is the analytical continuation of the series 



CO 

Tjk a k V {i+k+i) xk absolutely convergent in the unit circle, ff, being the coefficients of 
the original series. Parameter A is choosen from the condition of the most rapid 
convergence of series (|l|), that is from minimizing the quantity |1 — ^^7~^ \, where 
L is the step of truncation and F L (g;a,b) is the L-partial sum for F(g;a,b). In 
the regular scheme 0, [7], || parameters a and b are related to the exact asymptotic 
values a and b . Since, in practice we deal with a piece of the series only, where 
the asymptotic regime might not be established, we vary parameters a and b in 
a neighbourhood of their exact values Our principle observation is that the 



result of processing Fi(g; a, b) exhibit very weak dependence on the transformation 
parameters a and b varying in a wide range. The dependence becomes weaker with 
the growth of the approximation order and the smaller is the parameter of expansion 
g the better this property holds. We put the stability of the result of processing with 
respect to variation of a and b into the foundation of our technique to summation of 
divergent series. Such an approach allows us to apply the transformation ([]]) even 
if the exact asymptotic behaviour of the series being processed is unknown. 

The detailed analysis of the formulated resummation scheme applyed to the 
model functions 



H9) 



~ 9xA dx 



k=0 



- 1 ) kl 2 V , 



~'(xd a 



1 



-dx 



fc=0 



-l) k k\k bo g k 



1 + gx 

was performed in Ref. [22]. The convergence of the process for the function S(g), 
g = 1 is presented in Fig.l. Let us demonstrate our summation approach by esti- 
mating the ground state energy E(g) for several simple quantum mechanical sys- 
tems. Consider first the isotropic anharmonic oscillator p3| with the Hamiltonian 
H = x 2 + gx A . We observe the same stability of the result of processing E(g) with 
respect to a and b as for the model functions. The ground state energy estimates 
for (7 = 1 depending on the approximation order are listed in Table I. 



TABLE I. Numerical estimates for the isotropic anharmonic oscillator ground 
state energy at g = 1. 



L 


8 


9 


10 


11 


12 


Exact value 


S(l) 


1.392376 


1.392357 


1.392344 


1.392349 


1.392351 


1.392352 



For L = 8 our estimate is closer by one order to the exact value than the number 
1.391655 ± 0.004562 found in Ref. [[24] on the basis of Wynn's e-algorithm. 

The ground state energy E(g) of the cubic anharmonic oscillator with the Hamil- 
tonian 

-(x +y J 



H 



+ -\x 
4 



+ 2(1 - 5)xY + y 4 ] 
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depending of the anisotropy parameter 5 and the value of the coupling constant g 
was estimated in Ref.||18||. Those calculations were based on the knowlege of the 
exact values of the asymptotic parameters. The ground state energy estimated on 
the basis of our approach proved to be very close to the values of Ref . |fl8[| . The 
results given by two different methods are listed in Table II. The stability of the 
result of processing E(g) (g/4 = 0.1, 5 = —2.5) with respect to the variation of the 
parameters a and b is shown in Fig. 2. 

We have also studied the ground state energy E(g) of the Yukawa potential 

Vg( x ) = ~-exj)(-gx). 

x 

The dependence E(g) presented in Fig.3 is in agreement with the exact results 
[24] , p5|l . The exact critical value of g when the bound state disappears is g c = 
1.190612.... Using Winn's e-algorithm [[24] gives g c = 1.1836. The best estimate 
for g c in the frame of our approach yields g c — 1.191. In Fig. 4 we demonstrate the 
behavior of curves E(b) for g = 1.11 close to g c where the domain of stability of 
the result of processing begins to dissipate (c.f. with Fig.2 for the cubic anharmonic 
oscillator where the parameter of expansion g is small). 

In the context of this report, it is interesting to study the convergence of the RG 
series of certain conformal feild models, for which the exact results are known, using 
our summation procedure. In two dimentions, summation of e-series is a difficult 
problem because the parameter of expantion is large (e = 2). However, application 
of our approach gives, as a whole, relatively good estimates. So, using five-loop e- 
expansions for 0(iV)-symmetric model [2£| and setting N = 2cos(^), for the Ising 



model (m = 3) we obtain v = 0.925, r] = 0.220, 7 = 1.650 while the exact theory 
predicts v = 1, 77 = 0.25, 7 = 1.75. Obviously, the difference between our results and 
the exact values does not exceed 8%. For the model describing polymers (m = 2), 
our estimates are even better: v = 0.747 (exact 0.75), rj = 0.190 (exact 0.208), and 
7 = 1.352 (exact 1.350). Above numbers are also in an accordance with the results 
of Refs. m where a substantially different method was employed. 

The analysis of the simple models fulfilled above enables one to apply the sum- 
mation approach introduced to finding numerical estimates of important physical 
quantities in a number of real models. For example, it was used for calculation of 
critical exponents for the iV-vector field models describing magnetic phase transi- 
tions in cubic and tetragonal crystals |2^, pTfl . It can be expected that the developed 
technique may be useful in such QCD and QED as well. 

We thank Professor D. I. Kazakov for interesting and useful discussions of the 
problem of summation of the field-theoretical series during the XI International 
Conference "PROBLEMS OF QUANTUM FIELD THEORY". We are grateful to 
Professor J. Zinn- Justin for important remarks concerning the reference list of the 
paper. 
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TABLE II. Estimates of the ground state energy of the cubic anharmonic oscil- 
lator for various anisotropy parameters 5, coupling constant g, and approximation 
order L. 



Results by Kleinert et al, g/4 = 0.1 



L\S 


-2,5 


-1,5 


-0,5 


0,5 


1,5 


7 


1,217107 


1,192033 


1,164803 


1,134735 


1,100604 


9 


1,217107 


1,192034 


1,164810 


1,134736 


1,100604 


11 


1,217107 


1,192035 


1,164810 


1,134739 


1,100604 


Our estimates, g/A = 0.1 


L\S 


-2,5 


-1,5 


-0,5 


0,5 


1,5 


12 


1,21705 


1,19203 


1,16480 


1,134730 


1,00600 



relative error < 0.001% 

relative deviation from the results of Ref.|T8| < 0. 
0<6< 60 0.5 <a<1.5 

Results by Kleinert et al, g/A = 1.0 



001% 



L\S 


-2,5 


-1,5 


-0,5 


0,5 


1,5 


7 


1,941172 


1,862806 


1,773888 


1,669172 


1,535454 


9 


1,941172 


1,862815 


1,773909 


1,669188 


1,535425 


11 


1,941180 


1,862823 


1,773924 


1,669199 


1,535418 


Our estimates, g/A = 1.0 


L\S 


-2,5 


-1,5 


-0,5 


0,5 


1,5 


12 


1,9411 


1,8627 


1,7731 


1,6691 


1,5363 



relative error < 0.06% 

relative deviation from the results of Ref. |18| < 
0<6< 60 0.5 <o<1.5 



05% 
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Fig. 1 : Convergence of the estimates for the function £ (g), g=1 . 

Fig. 2: Resulting curves for the ground state energy E(g) of the cubic 
anharmonic oscillator from the 12-th approximation order. 
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Fig.3 
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Fig.4 
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Fig.3: Dependence of the ground state energy E(g) for the Yukawa potential 
from the 8-th approximation order. 



Fig.4: Resulting curves for the Yukawa potential ground state energy for g=1 .1 1 . 



